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I. INTRODUCTION 


The e“'"e“ ^ pp cross-section and the proton form fac¬ 
tor can be measured over a range of center-of-mass en¬ 
ergies by studying the initial state radiation (ISR) pro¬ 
cess e+e" —> pp'y (Fig. 1). The emission of a photon in 
the initial state gives rise to the possibility of measuring 
the cross-section of the nonradiative process e“'"e“ —> pp 
over a range of effective center-of-mass energies, from the 
threshold m = 2mp = 1.88 GeV/c^ to the full e+e” 
center-of-mass energy {^/s). The Born cross-section for 
this process, integrated over the nucleon momenta, is 
given by 




p—- 


dm d cos 9* 


= ‘^W{s,x,9*)app{m), ( 1 ) 


where m is the pp invariant mass, x = 2E*l^/s = 
1 — m?I and E* and 9* are the ISR photon energy 
and polar angle, respectively, in the e^e“ center-of-mass 
frame^. The function lT(s,a;,0*) [1], 


W{s,x,9;) 


a I 2 —2x + x"^ x^ 
TTx I sin^0* 2 


( 2 ) 


is the probability of ISR photon emission for 0* ^ 
m^j \fs^ where a is the fine-structure constant and me is 
the electron mass. The cross-section for the e+e” —*■ pp 
process is given by 


(Jpp{m) 


3m^ 


\GM{m)\'^ 


2^ 

m? 


|G£;(m)p 


( 3 ) 


with P = Y^l — Am^jm^, C = y/{l — e ^), and y = 
Tramp/{Pm) is the Coulomb correction factor [2], which 
makes the cross-section nonzero at threshold. The cross- 
section depends on the magnetic form factor (Gm) and 
the electric form factor (G_b); at threshold, \Ge\ = \Gm\- 
The modulus of the ratio of electric and magnetic form 
factors can be determined from the distribution of 0p, the 
angle between the proton momentum in the pp rest frame 
and the momentum of the pp system in the e“'"e“ center- 
of-mass frame. This distribution can be expressed as a 
sum of terms proportional to |GmP and \Ge\^- The full 


*Also with the Johns Hopkins University, Baltimore, Maryland 
21218 , USA 

^Also at Laboratoire de Physique Corpusculaire, Clermont- 
Ferrand, France 

^ Also with Universita di Perugia, Dipartimento di Fisica, Perugia, 
Italy 

^Also with Universita della Basilicata, Potenza, Italy 
^Deceased 

^ Throughout this paper, the asterisk denotes quantities in the 
e'*‘e“ center-of-mass frame. All other variables except 6p and 
are defined in the laboratory frame. 



differential cross-section for e“''e ^ ppy can be found, 

for example, in Ref. [3]. The 9p dependencies of the 
Ge and Gm terms are reminiscent of the sin^ 9p and 
1 -f cos^ 9p angular distributions for electric and magnetic 
form factors in the e“''e“ —*■ pp process. 

Measurements of the e“'"e“ —> pp cross-section have 
been performed in e“'"e“ experiments [4-9] with (20- 
30)% precision. The cross-section and proton form factor 
were deduced assuming |G£:| = |Gm|, and the measured 
proton angular distributions [5, 6] did not contradict this 
assumption. More precise measurements of the proton 
form factor have been performed in pp —> e“'"e“ experi¬ 
ments [10-12]. In the PSI70 experiment [10] at LEAR, 
the proton form factor was measured from threshold {pp 
annihilation at rest) up to a mass of 2.05 GeV/c^. The 
ratio \Ge/Gm\ was measured using the angular depen¬ 
dence of the cross-section and was found to be compatible 
with unity. The LEAR data show a strong dependence of 
the form factor on pp mass near threshold, and very lit¬ 
tle dependence in the range 1.95-2.05 GeV/c^. Analyses 
from Fermilab experiments E760 [II] and E835 [12] show 
a strong decrease in the form factor at higher masses, 
in agreement with perturbative QCD, which predicts a 
op,{rr?)lmP dependence. 

This work is an independent measurement by the 
BABAR Collaboration of the e+e" ^ pp cross-section 
(jpp{m), for pp masses up to 4.5 GeV/c^, based on the 
ISR process in e“*"e“ annihilation at a hxed center-of- 
mass energy near 10.6 GeV. This study significantly im¬ 
proves the measurement of app{m) in the pp mass range 
up to 3 GeV/c^. In contrast to previous e+e” and pp 
experiments, our measurement does not use the assump¬ 
tion that 1 G_b1 = [GmI- The ISR approach provides full 
9p coverage and hence high sensitivity to 1 G£;/Gm|- In 
this work, the mass dependence of the form-factor ratio 
\Ge/Gm\ is measured for pp masses below 3 GeV/c^. We 
also study J/p and p>{2S) production in e“''e“ ^ ppy, 
and measure the products F('0 —> e'^e~)B{'tp pp). 
A search for production of the y(4260) resonance, re¬ 
cently observed by BABAR in the ISR process e+e” —> 
E(4260)7 ^ J/'07r+7r“7 [13], is performed. 
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II. THE BABAR DETECTOR AND DATA 
SAMPLES 

We analyse a data sample corresponding to 232 fb“^ 
recorded with the BABAR detector [14] at the PEP-II 
asymmetric-energy storage ring. At PEP-II, 9-GeV elec¬ 
trons collide with 3.1-GeV positrons at a center-of-mass 
energy of 10.6 GeV (the r(4S) resonance). 

Charged-particle tracking is provided by a five-layer 
silicon vertex tracker (SVT) and a 40-layer drift chamber 
(DCH), operating in a 1.5-T axial magnetic field. The 
transverse momentum resolution is 0.47% at 1 GeV/c. 
Energies of photons and electrons are measured with a 
CsI(Tl) electromagnetic calorimeter (EMC) with a reso¬ 
lution of 3% at 1 GeV. Charged-particle identification is 
provided by specific ionization {dE / dx) measurements in 
the SVT and DCH, and by an internally reflecting ring¬ 
imaging Cherenkov detector (DIRC). Muons are identi¬ 
fied in the solenoid’s instrumented flux return, which con¬ 
sists of iron plates interleaved with resistive plate cham¬ 
bers. 

Signal and background ISR processes are simulated 
with Monte Carlo (MG) event generators based on 
Ref. [15], with the differential cross-section for e+e“ —> 
PP 7 taken from Ref. [3]. Because the polar-angle distri¬ 
bution of the ISR photon is peaked near 0° and 180°, 
the MG events are generated with a restriction on the 
photon polar angle: 20° < 9* < 160°, where 9* is mea¬ 
sured in the e“'"e“ center-of-mass frame. The extra soft- 
photon radiation from the initial state is generated with 
the structure function method [16]. To restrict the max¬ 
imum energy of the extra photons, the invariant mass 
of the hadron system combined with the ISR photon 
is required to be at least 8 GeV/c^. For background 
e“'"e“ —> 7 r“'' 7 r“ 7 , and processes, final 

state Bremsstrahlung is generated using the PHOTOS 
package [17]. Background from e“'"e“ —> qq is simulated 
with the JETSET [18] event generator. The response of 
the BABAR detector is simulated using the GEANT4 [19] 
program. The simulation takes into account the variation 
of the detector and accelerator conditions, and beam- 
induced background photons and charged particles over¬ 
lapping events of interest. 

III. EVENT SELECTION 

The preselection of e“'"e“ —> pp^ candidates requires 
that all the final-state particles are detected inside a fidu¬ 
cial volume. Since a significant fraction of the events 
contain beam-generated spurious tracks and photon can¬ 
didates, we select events with at least two tracks with 
opposite charge and at least one photon candidate with 
E* > 3 GeV. The polar angle of the photon is required 
to be in the well-understood region of the calorimeter: 
21.5° < 9^ < 137.5°. The charged tracks must origi¬ 
nate from the interaction point, have transverse momen¬ 
tum greater than 0.1 GeV/c, and be in the angular re¬ 


gion between 25.8° and 137.5°, so that particle identifica¬ 
tion (PID) may be performed using the DIRC detector. 
To suppress background from radiative Bhabha events, 
events in which each of the two highest momentum tracks 
has a ratio of calorimetric energy deposition to momen¬ 
tum in the range 0.9 to l.I are rejected. 

For events passing the preliminary selection, a kine¬ 
matic fit is performed to the e+e" —> hypoth¬ 

esis with requirements of total energy and momentum 
conservation. Here C can be e, p, tt, K or p, and 7 is the 
photon candidate with the highest energy in the e“'"e“ 
center-of-mass frame. For events with more than two 
charged tracks, the fit uses the two oppositely charged 
tracks that pass closest to the interaction point. The 
Monte Garlo simulation does not accurately reproduce 
the shape of the photon energy resolution function. This 
leads to a difference in the distributions of the of fh® 
kinematic fit for data and for MC simulated events. To 
reduce this difference, only the measured direction of the 
ISR photon is used in the fit; its energy is treated as a 
free fit parameter. For each of the five charged-particle 
mass hypotheses, the corrected angles and energies of the 
particles and the of the kinematic fit are calculated. 

The selection of e+e" —> ppj events relies upon both 
particle identification and event kinematics. The ex¬ 
pected number of events from the background processes 
e“'"e“ —> 7 r“'" 7 r“ 7 , pApL~x, and K~x significantly ex¬ 
ceeds the number of signal events (by two to three or¬ 
ders of magnitude). To suppress these backgrounds, both 
charged particles must be identified as protons accord¬ 
ing to the specific ionization (dE/dx) measured in the 
SVT and DGH, and the Gherenkov angle measured in the 
DIRG. These particle identification requirements lead to 
a loss of approximately 30% of the signal events, while 
suppressing backgrounds by factors of 15 x 10^, 500 x 10^, 
and 2 x 10 ^ for pion, muon, and kaon events, respectively. 

Background is further suppressed through require¬ 
ments on the of the kinematic fit: Xp < 30 and 
xl: > 30, where Xp and x^k are the x^ of the kinematic 
fit for the proton and kaon mass hypotheses, respectively. 
The distribution of Xp for Monte Carlo simulated ppj 
events is shown in Fig. 2 (left). The long tail in the dis¬ 
tribution at high x^ is due to events with extra photons 
emitted in the initial state. The dashed histogram is 
the Xp distribution for K^K~x Monte Carlo simulated 
events. Figure 2 (right) shows the distributions of xi- 
for K^K~ ■7 and pprf Monte Carlo simulated events with 
Xp < 30. The x^ requirements lead to a loss of 25% 
of signal events but provides additional background sup¬ 
pression by a factor of 50 for e“'"e“ ^ tt+tt^x and p~x 
events, and a factor of 30 for e“'"e“ —> K^K~x events. 

The pp invariant mass distribution is shown in Fig. 3 
for the « 4000 events that satisfy all the selection criteria. 
Most of the events have invariant mass below 3 GeV/c^. 
Glear signals from J/ijj ^ pp and 'tp{2S) pp decays are 
evident. 




tribution of the FSR amplitude is estimated to be [20] 
da/dm « |Faxp8mQ;^/3/(27s^), where Fax is the axial 
proton form factor. Assuming [Fax] Ri [Gm], the ratio 
of FSR to ISR cross-sections is determined to be about 
10“^ for pp masses below 4.5 GeV/c^, implying that the 
FSR background is sufficiently small to neglect. 


A. e“*“e — > 7r“*"7r 7, e"*“e —> K 7, 

e“*“e“ — > /i“*“/x“7 and e"*“e“ — > e“*'e“7 backgrounds 


FIG. 2: The Xp distribution (left) for MC simulated e^e“ —> 
ppX (solid line) and e'*"e“ ^ K^K~'y (dashed line) events, 
and the the distribution (right) for MC simulated e'’"e“ —> 
K^K~'y (solid line) and e'*'e“ —> pp7 (dashed line) events 
with Xp < 30. 



2 3 4 

Mpp(GeV/c") 


FIG. 3: The pp invariant mass spectrum for ppx candidates 
that satisfy all selection criteria. 


IV. BACKGROUND AND ITS SUBTRACTION 

The possible sources of background in the sample of 
e"'"e“ —> ppx candidates that pass the selection crite¬ 
ria described in the previous section include e“'“e“ —> 
7r+7r“7, e“'"e“ —> e'^e~ —> pd 'and e“'"e“ —*■ 

e"'"e “7 events in which the charged particles are misiden- 
tified as protons. Backgrounds from processes with pro¬ 
tons plus neutral particle(s) in the final state are also 
anticipated: e+e" —> pp 7 r°, pprj, ppir^j, etc. 

Of particular interest is the possible background from 
the process e“'"e“ —> pp') with the photon emitted from 
the final state. Due to different charge parity of the am¬ 
plitudes corresponding to initial state radiation and final 
state radiation (FSR), their interference does not con¬ 
tribute to the total e+e” ^ pp^ cross-section. The con¬ 


To estimate the background contribution from e“'“e“ ^ 
7 r+ 7 r“ 7 , data and Monte Carlo simulated events are se¬ 
lected with the following requirements on PID and on the 

of the kinematic fits: 

1 . one proton candidate, xj < 20 ; 

2. one proton candidate, Xp < 30, Xx > 30; 

3. two proton candidates, x^ < 20; 

4. two proton candidates, Xp < 30, x?k > 30. 

Here xj is the x^ of the kinematic fit for the pion mass 
hypothesis. 

The fourth set of conditions corresponds to the stan¬ 
dard selection criteria for ppx candidates. The invari¬ 
ant mass of the two charged particles under the 

pion-mass hypothesis is calculated; the Mn-ir distribu¬ 
tions for data selected with criteria 2 and 4 are shown 
in Fig. 4. The p resonance in the e“'"e“ —> 7 r“'' 7 r “7 re¬ 
action is clearly seen in the distribution corresponding 
to selection 2 (left plot in Fig. 4. The number of ttttx 
events with 0.5 < < 1 GeV/c^ passing each set of 

selection criteria is determined by fitting the distri¬ 
bution with a TTTry spectrum predicted by Monte Carlo 
plus a first order polynomial to account for background 
from non-TTTry processes. The Monte Carlo TTTry spectrum 
uses a model of the pion form factor based on existing 
experimental data. The results of the fits for 71717 can¬ 
didates passing selection criteria 1, 2, 3 and 4 are listed 
in Table I together with the corresponding numbers from 
the Monte Carlo simulation. 

Particle identification for the simulated 71717 events is 
accomplished using two sets of information: fully simu¬ 
lated observables that are used for particle identification 
in the same manner as in the analysis of data, and event 
weights for the simulated events based on pion misiden- 
tification rates derived from a control sample of known 
pions in data. The identification based on event weights 
does not take into account possible correlations between 
pion misidentification probabilities for two particles that 
overlap in the detector, or are in close proximity and 
therefore may underestimate the yield of wrongly identi¬ 
fied TTTry events. No events passing selection 4 are found 
in the fully simulated particle-ID sample, and so a 90% 
confidence level (CL) upper limit is estimated for the 
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TABLE I: The numbers of 7r7r7 events with 0.5 < < 1 GeV/c^ passing different selection criteria for data and Monte Carlo 

simulation (MC). WMC denotes Monte Carlo simulation with data-derived particle identification weights. The data numbers 
are obtained from the fit of the M-n-w distributions as described in the text. Rt^t: is the ratio of the numbers of events in the 
previous two rows. 



data 

1 proton candidate 
MG 

WMC 

data 

2 proton candidates 
MC 

WMC 

N{xi < 20 ) 

16200 ± 200 

21020 ± 230 

12300 ± 300 

190 ± 30 

246 ± 25 


35.5 ±0.8 

A'(Xp < 30,xx > 30) 

460 ± 120 

590 ± 40 

300 ±5 

- 

< 5.7 

0.90 ± 0.03 

Rtttt 

35 ±9 

36 ±2 

43 ± 1 

- 

> 43 


39 ±2 




FIG. 4: The spectrum for data events with Xp < 30 and 
Xrc > 30, and one proton candidate (left plot; selection 2 in 
the text) or two proton candidates (right plot; selection 4 in 
the text). The histograms are the results of the fit described 
in the text. 


standard selection (selection 4) in Table I. Because nei¬ 
ther the fully simulated nor weighted BID samples pre¬ 
dict the number of ttttx events seen in the data passing 
selection 4, an estimate is made based on the number of 
data events passing selection 3: N 4 = where 

Rtt-k is the ratio of the number of candidates that satisfy 
xi < 20 to the number that satisfy Xp < 30 and Xk > 30 
(the numbers given in the first and second rows in Ta¬ 
ble I. The statistical uncertainty on the scale factor Rt^^^ 
from the simulation is about 20%. Rt^t^ estimated for 
events with one and two misidentified pions are consis¬ 
tent with each other. Accordingly, the scale factor ratio 
R^^ = 35 ± 9 obtained from data is used, with an addi¬ 
tional 30% systematic uncertainty assigned. Finally, A ^4 
is estimated as N 4 = (190±30)/(35±14) = 5.4±2.3. The 
fit with A ^4 = 5.4, shown in Fig. 4, describes the mass dis¬ 
tribution for selection 4 very well. The total number of 
TTTry events remaining for the standard selection criteria 
is calculated as = 1.1 x N 4 = 5.9 ± 2.5, where 1.1 
is the ratio of the total number of ttttx events to those 
in the 0.5 < < 1 GeV/c^ mass region and is taken 

from simulation. The expected Mpp spectrum for TTTry 
events passing the ppj selection criteria is shown as the 
dotted histogram in Fig. 5. 

The procedure used to estimate the background from 
the e+e” —> K'^K~'y process is similar to that used to 
estimate the e“'"e“ ^ 7r“''7r“7 background. The number 
of events in the (j) meson peak in the distribution of in- 



2 2.2 2.4 2.6 

Mp-(GeV/c") 


FIG. 5: The calculated Mpp spectra for e'*'e“ ^ 7r’'"7r“7 
(dotted histogram), e''"e“ —> K~'y (solid histogram), and 
e^e“ —> /i^/r “7 (dashed histogram) background processes. 
The spectra are normalized to the number of events expected 
to pass the ppx selection criteria for each process: 5.9 ± 2.5 
for pions and 2.5 ± 1.0 for kaons. For the muon channel, the 
upper limit of 11 events is used for the normalization. 


variant mass of the charged particles calculated under 
the kaon hypothesis is used to determine the number of 
kaon events. The total number of KK^ events remain¬ 
ing after the standard selection criteria is estimated to 
be Nkk = 2.5 ± 1.0. The expected Mpp distribution for 
these events is shown as the solid histogram in Fig. 5. 

To estimate the electron background, the kinematic 
properties of the e“''e“ ^ e^e~x process are used. About 
60% of e“'"e “7 events have e“'"e“ invariant mass between 
3 and 7 GeV/c^ and costp* < —0.97, where tp* is the 
angle between the two tracks in the e“'"e“ center-of-mass 
frame. In the event sample with two proton candidates, 
only one event has the above characteristics. With this 
event assumed to be background from e“'"e“ ^ e“'"e“ 7 , 
the total e“''e “7 background is estimated to be 1.8 ± 1.8 
(0.8 events with Mpp < 4.5 GeV/c^). 
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The method used to estimate background relies 

on the difference between the two-proton and the two- 
muon mass spectra. From the simulation 44% of ii~^ 
events are expected to have a two-proton invariant mass 
greater than 4.5 GeV. In data, only four such events are 
found, with an expected background of 5 ± 3 events from 
the e+e“ ^ ppir^ process (see Sec. IVB). From these 
numbers, the total muon background is estimated to not 
exceed 11 events. A similar limit is obtained directly from 
Monte Carlo simulation. From about 2 million 
simulated events (20% of the number of events 

expected in data), no events pass the ppj selection crite¬ 
ria, leading to a 90% CL upper limit of 12 events. The 
expected Mpp spectrum for e+e” —*■ /r+/r“7 events nor¬ 
malized to a total of 11 events is shown as the dashed 
histogram in Fig. 5. This upper limit on the number of 
muon events is used as a measure of the systematic un¬ 
certainty due to pL^ background. This uncertainty is 
calculated as a function of the pp mass and is added to 
the systematic error on the number of ppj events. 


B. e'’"e —> ppiv^ background 

A dominant source of background to the e+e“ —> pp^ 
process arises from e^e“ —> ppir^. A significant fraction 
of ppTT^ events with an undetected low-energy photon or 
with merged photons from the decay are reconstructed 
under the pp"f hypothesis with a low value of and thus 
are not easily separable from the process under study. 
Experimental data is used to devise a procedure to sub¬ 
tract this background. 

For the ppir'^ background study, events with two 
charged particles identified as protons and at least two 
photons with energy greater than 0.1 GeV, one of which 
must have center-of-mass energy above 3 GeV, are se¬ 
lected. The two-photon invariant mass is required to 
be in the range 0.07 to 0.2 GeV/c^. A kinematic 
fit under the e+e” —> PP77 hypothesis is then per¬ 
formed. For events with more than two photons, all two- 
photon combinations are analyzed and only the combi¬ 
nation with the smallest iii the kinematic fit is con¬ 
sidered. Requirements on the of the kinematic fit 
(x^ < 25) and the two-photon invariant mass (0.1025 < 
Myj < 0.1675 GeV/c^) are then imposed on the e+e“ ^ 
ppTT^ candidates. The sidebands 0.0700 < < 

0.1025 GeV/c2 and 0.1675 < < 0.2000 GeV/c^ are 

used to estimate background. The Mpp spectra and cos 0p 
distributions for data events from the signal and sideband 
regions are shown in Fig. 6. The total number of selected 
events is 74 in the signal region and 10 in the sidebands. 
The number of e+e“ —*■ ppK^ events in the sidebands 
expected from MG simulation is 2.7. 

The PP77 selection criteria described above are applied 
to simulated e+e“ ^ qq events generated with the JET- 
SET package. The predicted number of e+e“ —> ppir'^ 
events is 73 ± 7. These events have an enhancement in 
the Mpp distribution near pp threshold, similar to that 




Mpp (GeV/cb cos Gp 

FIG. 6: The Mpp spectrum (left) and the cos6p distribu¬ 
tion (right) for selected e^e“ —> ppn° candidates in data. 
The shaded histogram shows the background contribution es¬ 
timated from Mjry sidebands. 



FIG. 7: The expected Mpp spectrum for e'’'e“ —> ppix^ events 
selected with the standard pp^ criteria. The spectrum is ob¬ 
tained by scaling the data distribution shown in Fig. 6 by the 
factor K-x[(j{Mpp) described in the text. 


in data (Fig. 6), but the angular distribution is peaked 
at cos 9p = ±1 and is not consistent with the nearly flat 
distribution found in data. To study these events, sim¬ 
ulated e+e” ^ pp7r° events are generated according to 
three-body phase space with an additional weight pro¬ 
portional to {Mpp — 1.86)2 (to imitate the Mpp distri¬ 
bution observed in data). The resulting cos 9p distribu¬ 
tion is flat. With these simulated events, Ky^(j{Mpp) 
is calculated as the ratio of the Mpp distributions for 
events selected with the standard ppry criteria to those 
with the ppTT^ criteria as a function of Mpp. The value 
of the ratio Ky[(^{Mpp) varies between 3.7 near Mpp 
threshold to 2.0 at 5 GeV/c^. The expected Mpp spec¬ 
trum for e“'"e“ ^ ppir^ background passing the ppj se¬ 
lection criteria is shown in Fig. 7 and is evaluated as 
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FIG. 8: Left: the ppTv° mass spectrum for selected e'’"e“ ^ 
ppit^'y candidates with < 20. The hashed histogram shows 
the background contribution estimated from the sideband 
with 30 < < 50. Right: the p-n^ (p7r°) mass spectrum 

for e^e“ ^ pp'K^'y data events with pp7r° mass outside the 
J/'i/i peak. 


■^MC(^pp) ^ {dN/dMpp)data, where {dN/dMpp)data is 
the mass distribution for e“'"e“ —> ppK^ events obtained 
above (Fig. 6). In Table II, the number of selected 
e+e" —> pp"f candidates and the expected number of 
e"'"e“ —*■ ppTT^ background events for different pp mass 
ranges are given. (The pp mass ranges near the J/'0 
and tpi^S) resonances are excluded.) The background 
contribution grows from 5% near pp threshold to 40% 
at Mpp « 4 GeV/c^. All observed pp'y candidates with 
Mpp >4.5 GeV/c^ are consistent with ppn^ background. 

The JETSET simulation is used to find other possible 
sources of background from e“''e“ —> qq. The number 
of qq events with hnal states other than ppn^ passing 
all cuts is 26 ± 4, with two final states, pp2'K^ and ppr;, 
accounting for 17 and 5 events, respectively. The back¬ 
ground contribution from these sources is estimated from 
data using the sideband as described below. 


C. e'’"e ^ pp7r°7 background 

The dominant ISR background process with protons 
in the final state is e+e” —> ppir^j. To estimate this 
background, events are selected with two charged parti¬ 
cles identified as protons and at least three photons with 
energy greater than 0.1 GeV, with one of these photons 
having center-of-mass energy above 3 GeV. The invari¬ 
ant mass of the two least energetic photons is required to 
be in the range 0.07-0.20 GeV/c^. For events that pass 
these criteria, a kinematic fit under the e“'"e“ ^ ppi:^^ 
hypothesis is performed. The distribution of ppi:^ invari¬ 
ant mass for events with x^ < 20 is shown in Fig. 8 (left). 
The shaded histogram shows the background contribu¬ 
tion estimated from the x^ sideband: 30 < x^ < 50. 
Most e“'"e“ —*■ ppir'^^ events have a ppK^ mass near a 
peak at 2.6 GeV. The contribution of the J/ip ^ ppn^ 
decay is also seen. The mass spectrum for 

events with ppn^ mass away from the J/ij) resonance is 
shown in Fig. 8 (right). The mass and width of the peak 


dominating in this distribution agree with the parameters 
of the V(1440) state, suggesting that the main mecha¬ 
nism in the e“'"e“ —*■ ppi:^ reaction is a transition through 
iV(1440)p or iV(1440)p intermediate states. 

The number of e+e“ —> ppK^j events passing the 
PP7 selection is estimated using Monte Carlo simulation; 
e+e” ^ ppTT^'y events are generated in the iV(1440)p7 -I- 
iV(1440)p7 model with the Np form factor reproducing 
the experimental ppm^ mass distribution. In the simu¬ 
lation, the ratio of detection efficiencies for the ppj and 
pp7r°7 selection criteria is (1.5±0.2)%. From 847±31 se¬ 
lected e+e“ —> ppTT^j candidates (Fig. 8) the background 
contribution to the sample of e+e” —> ppry candidates is 
estimated to be 13 ± 3 events (about 0.3% of the total 
number of selected ppj candidates). 

The background contribution from ISR processes with 
higher multiplicity is significantly lower. A procedure 
similar to that described above is used to estimate the 
background from the e^e“ —> pp2TT^^ process. Perform¬ 
ing a kinematic fit under the e+e" —> pp27r°7 hypothesis, 
560 ± 30 events are selected. From the Monte Carlo sim¬ 
ulation, the ratio of detection efficiencies for the pphf and 
pp27r°7 selection criteria is (0.09 ± 0.06)%, and the back¬ 
ground contribution due to e“'"e“ ^ pp2Tr^^ is estimated 
to be 0.5 ± 0.3 events. 


D. Background subtraction 


Table III summarizes the expected number of back¬ 
ground events estimated in the above sections. The uds 
column shows the number of background events expected 
from e^e“ —*■ qq with the ppn^ final state excluded. This 
background is estimated using the JETSET event gener¬ 
ator. Because JETSET has not been precisely verified 
for the rare processes contributing to the ppj candidate 
sample, the background estimation is based on the differ¬ 
ence in x^ distributions for signal and background events. 
The second row in Table III lists /?,, the ratio of N 2 , 
the number of events with 30 < Xp < 60, to Ni, the 
number of events with Xp < 30, calculated for signal 
and background processes using the Monte Garlo sim¬ 
ulation. The last row in Table III shows the expected 
numbers of signal and background events in the x^ side¬ 
band (30 < Xp < 60) calculated as N 2 = PNi. In the 
Table, it is evident that x^ distributions for signal events 
and those for background from the processes with higher 
hadron multiplicity (columns labeled uds and ppn^j) are 
very different. This difference can be used to estimate 
the background from these two sources, as follows. First, 
the 7r+7r“7, AT+AT^y, e+e^y, and ppn^ background de¬ 
termined in previous sections is subtracted from data. 
Then, from the resulting numbers of events in the signal 
and sideband x^ regions, N[ and N 2 , the numbers of sig¬ 
nal and background (from uds and ppir^j sources) events 
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TABLE II: The number of selected ppy candidates, Nppj, and the number of background events from the e^e —» ppn^ process, 
Npp^o, for different ranges of Mpp. The pp mass ranges near the J/ij) and V’(2<S') resonances are excluded. 


Mpp (GeV/c=) 

< 2.50 

2.50-3.05 

3.15-3.60 

3.75-4.50 

> 4.5 

^pp'y 

3166 

322 

37 

20 

4 

N - 0 

171 ± 29 

33± 11 

17 ±7 

8±4 

5±3 


TABLE III: Ni and N 2 are the numbers of selected pp'y candidates with a kinematic fit Xp < 30 and 30 < Xp < 60, respectively, 
for signal and for different background processes. The last column shows the numbers of candidates selected in data, fii is the 
ratio N 2 /N 1 obtained from simulation. The numbers for e''"e“ —> ppy are obtained from data using the background subtraction 
procedure described in the text. 



TT + TT y 

K+R-j 

ppyr® 

pp7r°7 

uds 

PP7 

data 

Ni 

5.9 ±2.5 

2.5 ± 1.0 

229 ± 32 

13 ±3 

26 ±4 

3737 ± 75 

4025 

f3i 

0.71 ±0.05 

0.52 ± 0.04 

0.13 ±0.01 

1.53 ±0.25 

1.44 ±0.30 

0.048 ± 0.001 


N 2 

4.2 ± 1.8 

1.3 ±0.5 

29 ±5 

20 ±3 

37 ±5 

179 ±5 

288 


with Xp < 30 can be calculated: 


corresponding bins (see Sec. VIII). 


N ■ — 


Ni - 


Npkn = N' - N,; 


^ — R - ! R. ’ “ -‘'I 

-L PvVl! Pbkg 


(4) 


V. ANGULAR DISTRIBUTION 


where Ptkg is the ratio of fractions of events in the side¬ 
band and signal x^ regions averaged over uds and 
backgrounds. For this coefficient (3bkg = 1.5±0.4 is used; 
it is the average of j3uds and with an uncertainty 

covering the full range of j3uds and fdppT^o^ variations. 

In Table III, it is also evident that pp^ events dominate 
the sideband. Therefore, the background is very sensitive 
to the accuracy of the /3pp.y coefficient. In particular, the 
data-Monte Carlo difference in the x^ distribution can 
lead to a systematic shift of the result. The simulation of 
the x^ distribution for ppy events is validated using data 
and simulated events in the channels e“'"e“ ^ pT^ and 
e+e" ^ K+K- '7, both of which are kinematically very 
similar to the process under study. In the simulations, 
the (3 coefficients for all three processes agree within 2%. 
The ratio the (3 coefficients for data and simulation is 
1.01 ± 0.03 for e+e" ^ K+R-j and 1.015 ± 0.012 for 
e+e" ^ pL~'). The ratio is used to correct 

the /3pp-Y value obtained from simulation, which results 
in (3pp^ = 0.048 ± 0.003. The error is estimated using the 
I3pp^ variation as a function of pp mass. 

With the method described above, the total number of 
e+e" —> PP7 events (Nsig) and background events from 
uds and ppn^j sources (Nbkg) in the signal region are 
3737 ± 67 ± 34 and 50 ± 12 ± 16, respectively. The main 
source of the systematic uncertainty on N^ig is the un¬ 
certainty in the ppK^ background. The numbers of uds 
and ppir^j background events are in good agreement with 
their estimations from simulation, (13 ± 3) -I- (26 ± 4) = 
39 ± 5. The total background in the Xp < 30 region is 
288 events, about 8% of the number of signal events. 

The background subtraction procedure is performed in 
each pp mass bin. The resulting numbers of signal events 
for each bin are listed in Table VI. The events from J/'0 
and '0(25') decays are subtracted from the contents of the 


The ratio of electric and magnetic form factors is ex¬ 
tracted by analysing the distribution of 9p, the angle be¬ 
tween the proton momentum in the pp rest frame and the 
momentum of the pp system in the e“'"e“ center-of-mass 
frame. In general, this distribution is given by 


dN _ 
d cos 9p 

-d ( Hm{ cos 9p, Mpp) -\- 


Ge 

Gm 


He{cos9p,Mpp) . (5) 


The functions Hm {cos 9p, Mpp) and He{cos9p, Mpp) do 
not have simple analytic forms, but are determined using 
MC simulation. Two samples of e“'"e“ —*■ ppj events are 
generated, one with Ge = d and the other with Gm = 
0. The obtained functions are similar to the 1 -I- cos^ 9p 
and sin^ 9p functions describing angular distributions for 
magnetic and electric form factors in the case of e+e” —> 
pp process. 

The angular distributions of the data are fit in six 
ranges of pp invariant mass from threshold to 3 GeV/c^. 
to measure \Ge/Gm\- The fit intervals, the correspond¬ 
ing numbers of selected events, and the estimated num¬ 
bers of background events are listed in Table IV. For each 
pp mass interval and each angular bin the background is 
subtracted using the procedure described in Section IV D. 
The angular distributions obtained are shown in Fig. 9. 
The distributions are fit to Eq. 5 with two free param¬ 
eters A (the overall normalization) and \Ge/Gm\- The 
functions Hm and He are modeled with the histograms 
obtained from MC simulation with the pp7 selection ap¬ 
plied. To account for differences between the pp mass dis¬ 
tributions of pprf events in data and MC simulation, the 
histograms Hm and He are re-calculated using weighted 
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TABLE IV: is the number of selected pp 7 candidates, Ni,kg 
is the number of background events, and \Ge/Gm\ is the 
fitted ratio of form factors, for each pp mass interval. 


Mpp, GeV/c^ 

N 

^bkg 

G_b/Gm| 

1.877-1.950 

533 

2 ± 7 

.n+0.24+0.17 
-*-•^-‘--0.22-0.12 

1.950-2.025 

584 

37± 12 

-1 

‘ O-0.25-0.14 

2.025-2.100 

602 

50± 15 

1 c:9+0-27+0.16 
1.9^-0,23-0.12 

2.100-2.200 

705 

42 ± 14 

1 1 q+0.20+0.12 

l.lO_o.i9-o,ll 

2.200-2.400 

592 

61 ± 16 

1 09+0.26+0.17 
-*--o-^_0.23-0.14 

2.400-3.000 

464 

45 ± 12 

1 99+0.30+0.16 
l-^^-0.30-0.16 








FIG. 9; The cos 6p distributions for different pp mass re¬ 
gions: (a) 1.877-1.950 GeV/c^ (b) 1.950-2.025 GeV/c^ 
(c) 2.025-2.100 GeV/c^, (d) 2.100-2.200 GeV/c^ (e) 2.200- 
2.400 GeV/c^ (f) 2.400-3.000 GeV/c^ The points with error 
bars show data distributions after background subtraction. 
The histograms are fit results: the dashed histograms show 
the contributions corresponding to the magnetic form factor; 
the dash-dotted histograms show the contributions from the 
electric form factor. 


events. The weights are obtained from the ratio of the pp 
mass distributions in data and simulation. In principle, 
the weights for Hm and He differ due to the different 
mass dependences of Gm and Ge- A first approxima¬ 
tion uses Gm = Ge- The fitted values of \Ge/Gm\ are 
then used in the next approximation to recalculate Hm 
and He- The second iteration leads to a small change 
(less than 2%) of the fitted values, and the procedure 
converges after a third iteration. 



-1 -0.5 0 0.5 1 

COS 0p 


FIG. 10: The angular dependence of the detection effi¬ 
ciency for simulated events with Mpp < 2.5 GeV/c^ before 
(open squares) and after (filled circles) correction for data- 
simulation difference in detector response. 

The simulated angular distributions are corrected to 
account for the differences between the data and the 
simulations, in particle identification, tracking, and pho¬ 
ton efficiencies. These corrections are discussed in de¬ 
tail in the next section. The angular dependences of 
detection efficiencies calculated with MC simulation be¬ 
fore and after corrections are shown in Fig. 10. The 
variations from uniform, which do not exceed 10% frac¬ 
tionally, derive from the momentum dependences of pro¬ 
ton/antiproton particle identification efficiencies. These 
manifest themselves as angular variations because there 
is a strong correlation between proton/antiproton mo¬ 
mentum and Op. In particular, the minima in detection 
efficiency at | cos0p| = 0.75 correspond to the minima in 
proton/antiproton identification efficiencies for momenta 
near 1.5 GeV/c. 

The histograms fit to the angular distributions are 
shown in Fig. 9; the values of \Ge/Gm\ are listed in 
Table IV and shown in Fig. 11. The curve in Fig. 11 
(l-|-aa:/(l-|- bx'^)) is used in the iteration procedure to 
calculate the weight. The quoted errors on \Ge/Gm\ 
are statistical and systematic. The dominant contribu¬ 
tion to the systematic error comes from the uncertainty 
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FIG. 11: The measured \Ge/Gm\ mass dependence. Filled 
circles depict BABAR data, the curve is the fit result. Open 
circles show the data from PS170 [10]. 


in the ppjr^ background. For example, for the 1.950- 
2.025 GeV/c^ range this contribution to the lower (up¬ 
per) error is 0.12 (0.16), which dominates the total sys¬ 
tematic error in this bin, and likewise dominates system¬ 
atic error in all bins. The error due to the limited MC 
simulation statistics (0.08 for 1.950-2.025 GeV/c^), un¬ 
certainties in the coefficients /3 used for background sub¬ 
traction (0.01), the uncertainty of description of mass de¬ 
pendence of |G_e/Gm| (0.01), and the uncertainty in the 
efficiency correction (0.02) are all considered. The last is 
conservatively estimated as the difference between fitted 
values of \Ge/Gm\ obtained with and without applying 
the efficiency correction. 

The angular distribution for J/ip ^ pp decay has also 
been studied. Its shape is commonly parameterized us¬ 
ing the form 1 -|- a cost'd. The coefficient a has been 
measured with relatively high precision in several experi¬ 
ments [24, 25], and its average value is a = 0.660±0.045. 
The BABAR data distribution for J/ip pp decay is 
shown in Fig. 12. The non-peaking background is sub¬ 
tracted by taking the difference between the histograms 
for the signal mass region (3.05-3.15 GeV/c^) and the 
mass sidebands (3.00-3.05 and 3.15-3.20 GeV/c^). The 
fitting procedure used is similar to the one described 
above with a = (1 — gr)/(l ± gr), where g = \Ge/Gm\^ 
and T = The resulting value a = O.TSto^as is 

in agreement with the world average but has significantly 
larger uncertainty. 

To cross-check this method to measure IG^/GmI, a 
comparison is made between the data and simulated dis¬ 
tributions of cosOk for e“''e“ (p^ ^ process. 

Here 9k is defined analogous to the definition of 6 p: 9 k is 



FIG. 12: The cos 9p distributions for J/ip pp decay. 
The points with error bars correspond to the background- 
subtracted data distribution, the solid histogram is the fit re¬ 
sult, and the dashed and dashed-dotted histograms show the 
fit contributions from the magnetic and electric form factors, 
respectively. 

the angle between the K~ momentum in the K'^K~ rest 
frame and the momentum of the K'^K~ system in the 
e“'"e“ center-of-mass frame. The angular dependence for 
this process is well known (approximately cos^ 9k) and 
event kinematics are similar to the e“'"e“ —> pp^ kine¬ 
matics near threshold. Figure 13 shows the ratio of data 
and simulated distributions over cos 9 k for events with 
the K+K- mass near the (p. The simulation describes 
the angular dependence of detection efficiency well. 

When these BABAR measurements of the |G£;/Gm| 
ratio are compared with the PS170 measurements [10] 
(Fig. 11), a large disagreement is seen for Mpp larger 
than 1.93 GeV/c^. 

VI. DETECTION EFFICIENCY 

The detection efficiency, determined using Monte Carlo 
simulation, is the ratio of true pp mass distributions com¬ 
puted after and before applying selection criteria. Be¬ 
cause the e“'"e“ ^ ppj differential cross-section depends 
on the form factors, the detection efficiency is somewhat 
model-dependent. The model used in this study has the 
\Ge/Gm\ ratio obtained from a fit of experimental angu¬ 
lar distributions (curve in Fig. 11) for Mpp < 3 GeV/c^, 
and |G_e/Gm| = 1 for higher masses. The detection ef¬ 
ficiency calculated in this model, shown in Fig. 14, is fit 
to a third-order polynomial for Mpp < 3 GeV/c^ and a 
constant for Mpp > 3 GeV/c^. The statistical error of 
the detection efficiency is about 1%. The model error 
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-1 -0.5 0 0.5 1 

COS 0j^ 


FIG. 13: The ratio of data and simulated distributions of 
cos 9 k for e'’"e“ ^ K'^K~'y process. 




3 3.5 4 4.5 

Mpp(GeV/c2) 


FIG. 14: The pp mass dependence of detection efficiency 
obtained from MC simulation. The line on the left plot 
is the fit to a third-order polynomial. The efficiency for 
Mpp > 3 GeV jc? is fit to a constant value. 


is determined from the uncertainty in the \Ge/Gm\ ra¬ 
tio: for Mpp < 3 GeV/c^, varying the ratio within its 
experimental uncertainty leads to a 1% change in the de¬ 
tection efficiency. This is taken as the model error. This 
small value is not surprising, due to the relatively small 
difference between the detector sensitivities for pure elec¬ 
tric and magnetic transitions. This difference was calcu¬ 
lated with simulated event samples in which Ge = 0 and 
Gm = 0 and is shown as a function of Mpp in Fig. 15. It 
does not exceed 20%. For masses above 3 GeV/c^, where 
the \Ge/Gm \ ratio is unknown, a 10% model error equal 
to half of the difference between detection efficiencies cor¬ 
responding Gm = 0 and Ge = 0 is used. 

The efficiency determined from MG simulation (emc) 
must be corrected to account for data-MG simulation dif- 
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FIG. 15: The relative difference between detection efficiencies 
for a purely electric (Gm = 0) and purely magnetic {Ge = 0) 
transition in e^e“ —> pp'y reaction. 


ferences in detector response: 

£ = £mC J^(l-b <^i)j (6) 

where 5i are efficiency corrections for each of several ef¬ 
fects. These corrections are discussed in detail below and 
summarized in Table V. 

Inaccuracy in the simulation of angular and momen¬ 
tum resolutions and radiative corrections may account for 
some of the data-MG difference in the fraction of events 
rejected by the requirement that Xp < 30. The efficiency 
correction for this effect is estimated by comparing data 
and simulated distributions for the e“'"e“ —> 
process, which has kinematics similar to the process un¬ 
der study. An exclusive e+e” ^ sample is se¬ 

lected by requiring that both charged tracks be identi¬ 
fied as muons. To remove possible background contri¬ 
butions from hadronic events with J/'0 —> p.~ de¬ 

cay, events with di-muon invariant mass in the range 
3.0 < Mpp < 3.2 GeV/c^ are excluded. The ratio 
of the number of selected muon events with > 30 
and Xp < 30 varies from 0.35 to 0.4 in the Mpp range 
from threshold to 4.5 GeV/c^. When comparing data 
and MC simulation in the region Xp > 30, the cut 
Mppj < 8 GeV/c^ needs to be applied, for consistency, 
to the data sample, since this cut is already applied into 
the MC simulation and therefore some events of this non¬ 
signal region are rejected in the simulated sample. To 
characterize data-MG simulation difference in the x^ dis¬ 
tribution, a double ratio (k) is calculated as the ratio 
of Af(Xp > 30)/Af(Xp < 30) obtained from data to the 
same quantity obtained from MC simulation. The value 
of the double ratio is k = 1.04 ± 0.01, essentially inde¬ 
pendent of mass. The efficiency correction for the x^ cut 
is calculated as 

jV(x^ < 30) + fV(x^ > 30) 

^ N{x^ < 30) -b kN{x^ >30) ’ ^ 

where N{x^ < 30) and N{x^ > 30) are the numbers of 
simulated ppx events with x^ < 30 and x^ > 30, re- 
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spectively. The values of the efficiency correction i5i for 
different pp masses are listed in Table V. Its statistical 
error is about 0.3%. An additional 1% systematic error, 
equal to the correction variation in the pp mass region of 
interest, is added in quadrature. 

The effect of the > 30 cut is studied using e+e“ ^ 
J/'f/’y —> PP7 events. The number of J/tp’s is determined 
using the sideband subtraction method. The event losses 
due to Xk > 30 cut are found to be 1.7 ± 0.7% in data 
and 1.7 ± 0.2% in MC simulation. As the data and sim¬ 
ulated values are in good agreement, there is no need to 
introduce any efficiency corrections for the > 30 cut. 
The systematic uncertainty associated with this cut is 
0.7%. 

Another possible source of data-MC simulation differ¬ 
ences is track loss. The systematic uncertainty due to dif¬ 
ferences in track reconstruction is estimated to be 1.3% 
per track. Specifically, for e+e" —> ppj only, the sys¬ 
tematic error can originate from slightly imperfect simu¬ 
lation of nuclear interactions of protons and antiprotons 
in the material before the SVT and DCH. The simula¬ 
tion shows that nuclear interaction leads to the loss of 
approximately 6% of e+e” —> ppx events. For data-MC 
simulation comparison, a specially selected event sample 
with A(A) decaying into p{p)n is used. The A are selected 
by imposing requirements on pir invariant mass and the 
A flight distance. The amount of material before the 
SVT (1.5% of nuclear interaction length) is comparable 
to the amount of material between the SVT and the DCH 
(1.4% of nuclear interaction length). The probability of 
track losses between the SVT and the DCH is measured 
from the A(A) sample. The data and simulation proba¬ 
bilities are found to be in good agreement for protons. A 
substantial difference is observed for antiprotons, which 
is consistent with a large (a factor of 2.6 ± 1.0) over¬ 
estimation of the antiproton annihilation cross-section in 
simulation. This difference in the antiproton annihilation 
cross-section in data and simulation leads to a correction 
of about (1.0 ± 0.4)% to the detection efficiency for ppry 
events. 

The data-MC simulation difference in the particle iden¬ 
tification is studied with use of events with a J/ip —^pp 
decay. Due to the narrow J/ip width and low back¬ 
ground, the number of J/'f/’ —*■ pp decays may be deter¬ 
mined using selections with either one or two identified 
protons. The background from non-J/ip events is sub¬ 
tracted using sidebands. The p/p identification probabil¬ 
ities are determined as functions of the p/p momenta by 
calculating the ratio of the number of events with both 
the proton and the antiproton identified to the number 
of events with only one identified proton or antiproton. 
The ratio of data-MC identification probabilities is used 
to reweight selected simulated events and calculate effi¬ 
ciency corrections. The correction is about (3±3)% and 
varies within 0.5% depending on pp mass. The error in 
the correction factor is determined from the statistical 
uncertainty in number of selected J/ip events. 

Another correction must be applied to the photon de¬ 


TABLE V: The values of different efficiency corrections 5i for 
pp masses 1.9, 3.0, and 4.5 GeV/c^. 


effect 


A(3),% 

<5i(4.5),% 

Xp < 30 cut 

Xk > 30 cut 

-0.7 ± 1.0 

-1.1 ± 1.0 

-1.7 ± 1.0 

0.0 ±0.7 

0.0 ±0.7 

0.0 ±0.7 

track reconstruction 

0.0 ±3.0 

0.0 ±3.0 

0.0 ±3.0 

nuclear interaction 

0.8 ±0.4 

1.1 ±0.4 

1.0 ±0.4 

PID 

2.5 ±3.3 

3.2 ± 2.4 

3.5 ±2.7 

photon inefficiency 

-1.3 ±0.1 

-1.3 ±0.1 

-1.3 ±0.1 

photon conversion 

0.4 ±0.2 

0.4 ±0.2 

0.4 ±0.2 

trigger 

-0.6 ±0.3 

- 

- 

total 

1.1 ±4.7 

2.3 ±4.1 

1.9 ±4.2 


tection efficiency. There are two main sources for this 
correction: data-MC simulation differences in the proba¬ 
bility of photon conversion in the detector material before 
the DCH, and the effect of dead calorimeter channels. A 
sample of e“'"e“ ^ pTx events is used to determine 
the photon inefficiency in data. Events with exactly two 
oppositely charged tracks identified as muons are selected 
and a kinematic fit is performed, constraining zero recoil 
mass against the muon pair. A tight cut on x^ of the kine¬ 
matic fit selects events with only one photon in the final 
state. The photon direction is determined from the fit. 
The photon detection inefficiency is calculated using the 
ratio of number of events not passing the > 3 GeV 
cut to the total number of selected p~x events. The 
obtained photon inefficiency 3.3% can be compared to the 
2% inefficiency in e“'“e“ —*■ ppy simulation. The observed 
data-MC difference in the photon inefficiency leads to an 
efficiency correction of (—1.3 ± 0.1)% that is practically 
independent of pp mass. The data-MC simulation differ¬ 
ence in the probability of photon conversion is studied 
using e+e“ ^ 77 events and found to be (0.4 ± 0.2)%. 

The quality of the simulation of the trigger efficiency 
is also studied. The overlap of the samples of events 
passing different trigger criteria and the independence of 
these triggers are used to measure the trigger efficiency. 
A small difference ((—0.6±0.3)%) in trigger efficiency be¬ 
tween data and MC simulation is observed for pp masses 
below 2.025 GeV/c^. 

All efficiency corrections are summarized in Table V. 
The corrected detection efficiencies are listed in Table VI. 
The uncertainty in detection efficiency includes a simu¬ 
lation statistical error, a model uncertainty, and the un¬ 
certainty of the efficiency correction. 


VII. e+e- ^ pp CROSS-SECTION AND 
PROTON FORM FACTOR 


The cross-section for e“'"e —> pp is calculated from the 

pp mass spectrum using expression 


(Jpp{m) 


(dN/dm) 

corr 

eRdL/dm ’ 


( 8 ) 
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where m is the pp invariant mass, {dN/dm)corr is the 
mass spectrum corrected for resolution effects, dL/dm 
is the ISR differential luminosity, e is the detection effi¬ 
ciency as a function of mass, and i? is a radiative correc¬ 
tion factor accounting for the Born mass spectrum dis¬ 
tortion due to emission of several photons by the initial 
electron or positron. The ISR luminosity is calculated us¬ 
ing the total integrated luminosity L and the probability 
density function for ISR photon emission in Eq. 2: 

dL a , 1 -I- C 9^\ 2m ^ , 

Here x = 1 —m^/s, -^/s is the e+e“ center-of-mass energy, 
C = cos 9q , and 9q determines the range of polar angles 
in the e+e“ center-of-mass frame: 9q < 9* < 180° — 9q 
for the ISR photon. In this study, 9q = 20°, because the 
detector efficiency is determined using a simulation with 
20° < 9* < 160°. The integrated ISR luminosity for each 
AIpp bin is listed in Table VI. 

The radiative correction factor R is determined us¬ 
ing Monte Carlo simulation at the generator level, with¬ 
out any detector simulation. Two pp mass spectra are 
generated; the first using the pure Born amplitude for 
the process e+e” —> ppj, and the second using a model 
with higher-order radiative corrections included with the 
structure function method [16]. The radiative correc¬ 
tion factor R is the ratio of the second spectrum to 
the first, and varies from 1.002 at the pp threshold to 
1.02 at 4.5 GeV mass. The value of R depends on 
the requirement on the invariant mass of the ppy sys¬ 
tem. The R in this study corresponds to the requirement 
Mpp-y > 8 GeV/c^ imposed in the simulation. The the¬ 
oretical uncertainty in the radiative correction calcula¬ 
tion with the structure function method does not exceed 
1% [16]. To check the theoretical uncertainty, a compar¬ 
ison of the cross-sections calculated with the structure 
function method and the Phokhara [3] event generator is 
performed. The Phokhara generator uses formulae with 
next-to-leading order radiative corrections in the initial 
state. The uncertainty of Phokhara generator is esti¬ 
mated to be less than 1% [26]. The ratio of pp mass 
spectra obtained with the two generators differs from 
unity by about 1% and does not contradict estimates 
of the theoretical uncertainties. The radiative correc¬ 
tions calculated include initial state radiation and the 
effect of loops at electron vertex, but do not include cor¬ 
rections for leptonic and hadronic vacuum polarization 
in the photon propagator. Cross-sections obtained with 
such corrections are sometimes referred to as “dressed” 
cross-sections while those which account fully for higher 
order processes are referred to as “bare” cross-sections. 
See Ref. [27] for a more complete discussion. 

The resolution-corrected mass spectrum is obtained by 
unfolding the mass resolution from the measured mass 
spectrum. Using the MC simulation, the migration ma¬ 
trix A is obtained, representing the probability that an 
event with true mass (MV“°) in bin j is reconstructed in 



FIG. 16: The e'’"e“ —> pp cross-section measured in this 
work and e^e“ experiments: FENICE[6], DM2[5], DM1[4], 
ADONE73[7], BES[8], CLEO[9]. The contribution of J/'ij} —> 
pp and '0(25) —> pp decays is subtracted. 



FIG. 17: The e"’"e pp cross-section near threshold 

measured in this work and e'*'e“ experiments: FENICE[6], 
DM2 [5], DM1 [4], ADONE73[7], BES[8]. 


bin i : 



As the chosen bin width significantly exceeds the mass 
resolution for all pp masses, the migration matrix is 
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TABLE VI: pp invariant mass {Mpp), number of selected events {N) with background subtracted, detection efficiency (e), ISR 
luminosity (L), measured cross-section (upp), and Fp, the effective form factor for e''"e“ —> pp. The contribution of J/tp —> pp 
and '!/>(2S) —> pp decays has been subtracted. The quoted uncertainties in N and a are statistical and systematic. For the form 
factor, the combined uncertainty is listed. 


Mpp (GeV/c^) 

N 


e 


L (pb-i) 

opp (pb) 

|Ap| 

1.877- 

-1.900 

157 ± 

13 ±3 

0.171 

± 

0.008 

1.141 

802 ± 68 ± 43 

0.453l°°i 

1.900- 

-1.925 

190 ± 

15 ±3 

0.173 

± 

0.008 

1.236 

887 ± 71 ± 46 

0.354l°;°l^ 

1.925- 

-1.950 

180 ± 

15 ±3 

0.175 

± 

0.008 

1.254 

819 ± 68 ± 43 

0.305l°:;)l® 

1.950- 

-1.975 

171 ± 

15 ±4 

0.177 

± 

0.008 

1.272 

760 ± 68 ± 41 

U.Z/D_0.015 

1.975- 

- 2.000 

176 ± 

16 ±5 

0.178 

± 

0.008 

1.290 

765 ± 68 ± 44 

U.ZOO_Q o;L5 

2 . 000 - 

-2.025 

201 ± 

17 ±5 

0.180 

± 

0.008 

1.308 

854 ± 71 ± 48 

n o'yo+O-OiS 
U.Z/O_0.014 

2.025- 

-2.050 

181 ± 

16 ±6 

0.182 

± 

0.008 

1.328 

748 ± 68 ± 45 

0.250t°“« 

2.050- 

-2.075 

196 ± 

17 ±6 

0.183 

± 

0.008 

1.346 

794 ± 68 ± 46 

0.254t°°« 

2.075- 

- 2.100 

180 ± 

17 ±5 

0.184 

± 

0.008 

1.365 

715 ± 66 ± 40 

U.Z,Oy_Q q;L 3 

2 . 100 - 

-2.125 

203 ± 

18 ±5 

0.185 

± 

0.008 

1.383 

794 ± 68 ± 44 

U.ZOU_0.013 

2.125- 

-2.150 

188 ± 

17 ±5 

0.185 

± 

0.009 

1.402 

721 ± 66 ± 40 

n 907 + 0.012 
U.ZrO/ _0.013 

2.150- 

-2.175 

147 ± 

15 ±5 

0.186 

± 

0.009 

1.421 

554 ± 58 ± 34 

n 007+0012 

U.ZU/_o.oi3 

2.175- 

- 2.200 

128 ± 

15 ±6 

0.187 

± 

0.009 

1.440 

477 ± 56 ± 31 

0.19lt°;°ll 

2 . 200 - 

-2.225 

119 ± 

14 ±5 

0.187 

± 

0.009 

1.459 

435 ± 52 ± 29 

0.183t°:;)12 

2.225- 

-2.250 

109 ± 

13 ±3 

0.187 

± 

0.009 

1.478 

392 ± 47 ± 21 

0 1 74 + 0-011 
U.i /+_o.oi 2 

2.250- 

-2.275 

69 ±: 

11 ±3 

0.188 

± 

0.008 

1.497 

245 ± 40 ± 16 

A 1 07 + 0.012 

U.J.O< _0.013 

2.275- 

-2.300 

70 ±: 

11 ±5 

0.188 

± 

0.008 

1.516 

244 ± 39 ± 20 

A 1 07 + 0.012 

U.iO/ _0.013 

2.300- 

-2.350 

82 ±: 

12 ±9 

0.188 

± 

0.008 

3.092 

140 ± 21 ± 16 

0 1 r)c;+0-009 

U.iUO_Q 0 ;LO 

2.350- 

-2.400 

80 ± : 

11 ± 7 

0.188 

± 

0.008 

3.172 

133 ± 19 ± 13 

n 1 OQ+O-OOS 
U.1UO_0.009 

2.400- 

-2.450 

91 ±: 

11 ±2 

0.187 

± 

0.008 

3.251 

149 ± 18 ± 7 

o.iiot°:“I 

2.450- 

-2.500 

52 ± 

9±2 

0.187 

± 

0.008 

3.331 

83 ± 15 ± 5 

0.083t°:““« 

2.500- 

-2.550 

63 ±: 

10 ±2 

0.186 

± 

0.008 

3.414 

100 ± 16 ± 6 

0 nQ2+0-007 

U.UyZ,_Q QQg 

2.550- 

-2.600 

39 ± 

8±2 

0.185 

± 

0.008 

3.496 

60 ± 13 ± 4 

0 072+0-007 
U.U/^_0.008 

2.600- 

-2.650 

31 ± 

8±2 

0.183 

± 

0.008 

3.580 

47 ± 11 ±4 

U.UOO_0.009 

2.650- 

-2.700 

25 ± 

7±2 

0.182 

± 

0.008 

3.664 

37± 10±4 

0 05Q+0-008 
U.UOy_Q QQQ 

2.700- 

-2.750 

21 ± 

7±2 

0.180 

± 

0.008 

3.749 

31 ±10 ±4 

0 nc;4+0-008 

U.UO+_0.010 

2.750- 

-2.800 

25 ± 

7±2 

0.179 

± 

0.008 

3.837 

37± 10±4 

U.UOU_Q o;LO 

2.800- 

-2.850 

20 ± 

6±2 

0.178 

± 

0.008 

3.924 

30 ± 9 ± 3 

U.UO+_0.009 

2.850- 

-2.900 

19 ± 

6±2 

0.176 

± 

0.008 

4.013 

27±9±3 

U.UOZ_Q 

2.900- 

-2.950 

19 ± 

6±2 

0.175 

± 

0.007 

4.103 

26 ± 8 ± 3 

0 0^2+0-008 
U.UOZ_Q 009 

2.950- 

-3.000 

9±! 

5±2 

0.173 

± 

0.007 

4.195 

12±7±3 

U.UOO_Q 014 

CO 

b 

O 

O 

-3.200 

11 ± 

9±8 

0.169 

± 

0.018 

17.719 

3.6 ±3.0 ±2.8 

0 021 +0-009 

U.UZi_0.021 

3.200- 

-3.400 

8±; 

5±7 

0.169 

± 

0.018 

19.289 

2.3 ± 1.6 ±2.0 

0 01 7+0-008 
U.Ul 1 _0.017 

3.400- 

-3.600 

6 ±' 

4±3 

0.169 

± 

0.018 

20.960 

1.7 ± 1.0 ±0.9 

0 01 0+0-005 
U.UiD_o.o09 

3.600- 

-3.800 

8 i ^ 

4±3 

0.168 

± 

0.018 

22.739 

2.2 ± 1.1 ± 0.8 

0 01 Q+0 005 

U.Uiy_0.008 

3.800- 

-4.000 

5±: 

3±3 

0.168 

± 

0.018 

24.645 

1.2 ±0.8 ±0.7 

0 01 c;+0-005 
U.UlO_0.009 

4.000- 

-4.250 

4±: 

3±3 

0.168 

± 

0.018 

33.701 

0.7 ±0.5 ±0.5 

0 . 0111 °;“?® 

4.250- 

-4.500 

i±: 

3±3 

0.167 

± 

0.018 

37.214 

0.1 ±0.4 ±0.5 

0.005l?;r5 

nearly diagonal, with the values of diagonal elements 
~ 0.9, and next-to-diagonal ~ 0.05. We unfold the mass 
spectrum by applying the inverse of the migration ma¬ 
trix to the measured spectrum. This procedure changes 

quoted are statistical and systematic. The systematic 
uncertainty includes the uncertainty in the number of 
signal events and detection efficiency, an error of total 
integrated luminosity (1%), and the uncertainty in the 


the shape of the mass distribution insignificantly, but in¬ 
creases the errors (by w20%) and their correlations. 

The number of events in each mass bin is listed in 
Table VI. The quoted errors are statistical and system¬ 
atic (with the systematic errors due to uncertainties in 
background subtraction). The calculated cross-section 
for e“''e“ —> pp is shown in Fig. 16 and listed in Ta¬ 
ble VI. For mass bins 3-3.2 GeV/c^ and 3.6-3.8 GeV/c^, 
the nonresonant cross-section is quoted with J/tp and 
ip{2S) contributions excluded (see Sec. VIII). The errors 


radiative corrections (1%). A comparison of this result 
with available e“'"e“ data is shown in Fig. 16 and the 
near-threshold region is shown in Fig. 17. 

The e“'"e“ —*■ pp cross-section is a function of two form 
factors, but due to poor determination of the \Ge/Gm\ 
ratio, they cannot be extracted from the data simultane¬ 
ously with reasonable accuracy. Therefore, the effective 
form factor is introduced: 



( 11 ) 
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FIG. 18: The proton effective form factor measured in this 
work and in e^e“ and pp experiments: FENICE[6], DM2[5], 
DM1[4], BES[8], CLEO[9], PS170[10], E835[12], E760[ll]. 
The upper plot shows the mass interval from pp threshold 
to 3.01 GeY/c?. The lower plot presents data for pp masses 
from 2.58 to 4.50 GeV/c^ 


FIG. 19: The proton effective form factor measured in this 
work and in e''"e“ and pp experiments, shown on a logarith¬ 
mic scale: FENICE[6], DM2[5], DM1[4], BES[8], CLEO[9], 
PS170[10], E835[12], E760[ll]. The curve corresponds to the 
QCD fit described in the text. 


where <jpp{m) is the measured e+e" —> pp cross-section 
and (Jnijn) is the cross-section obtained from Eq. 3 
under the assumption that \Ge\ = \Gm\ = 1- At 
Mpp = 2 GeV/c^ cr„ ~ 10 nb. This definition of 
the effective form factor Fp{m) permits comparison of 
our measurements with measurements from other exper¬ 
iments, in e“'"e“ as well as pp collisions. Most available 
form-factor data are analyzed using the assumption that 
\Ge\ = \Gm\- The calculated effective form factor is 
shown in Fig. 18 (linear scale), in Fig. 19 (logarithmic 
scale) and in Table VI. The form factors here are av¬ 
eraged over bin width, and the four points of PS170 [10] 
with lowest mass are all situated within the first bin of the 
BABAR measurement. For the mass region near thresh¬ 
old where the form factor changes rapidly with mass, 
the cross-section and effective form factor with a smaller 
bin size are calculated. These results are listed in Ta¬ 
ble VII. The effective form factor is shown in Fig. 20. 
In Figs. 18, 19 and 20, it is evident that the BABAR 
effective form factor results are in reasonable agreement 
with those of other experiments. The form factor has 
a complex mass dependence. The significant increase in 
form factor as the pp threshold is approached may be a 
manifestation of a pp subthreshold resonance [21]. The 
rapid decreases of the form factor and cross-section near 
2.25 GeV/c^ and 3 GeV/c^ have not been discussed in 
the literature. The dashed line in Fig. 19 corresponds 
to the asymptotic QCD fit [22] for proton form factor 
Fpp ~ ~ C/(m^ log (m^/A^)), applied to all 



Mpp(GeV/c") 


FIG. 20: The proton effective form factor near pp thresh¬ 
old measured in this work and in e^e“ and pp experiments: 
FENIGE[6], DM1[4], PS170[10]. 
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TABLE VII: pp invariant mass {Mpp), number of selected events {N) after background subtraction, measured cross-section 
(cTpp), and effective form factor for e^e“ —> pp. The quoted errors in N and app are statistical and systematic. For the effective 
form factor, the total combined error is listed. 


Mpp (GeV/c=) 

N 

app (pb) 

I^pI 

1.8760-1.8800 

18 ± 5 ± 1 

656 ± 161 ± 40 

0.574l°;°S 

1.8800-1.8850 

34 ± 6 ± 1 

808 ± 155 ± 43 

u.^yo_Q Q52 

1.8850-1.8900 

27±6± 1 

656 ± 154 ± 36 

u.oyu_o.o50 

1.8900-1.8950 

37±7±1 

889 ± 174 ± 48 

0.419l°;°S 

1.8950-1.9000 

38 ± 8 ± 1 

901 ± 182 ± 48 

0.3980^44 

1.9000-1.9050 

42 ± 9 ± 1 

995 ± 207 ± 56 

U.Oyy_Q Q4g 

1.9050-1.9100 

31 ±8± 1 

726 ± 186 ± 41 

u.ozu_Q Q4e 

1.9100-1.9150 

49 ± 9 ± 1 

1138 ±210 ±60 

D QQ7+0-036 
u.oy < -0.040 

1.9150-1.9250 

69 ± 10 ± 1 

798 ± 116 ±43 

n *^91 +0-024 

U.OZl_Q Q26 

1.9250-1.9375 

91 ± 11 ± 2 

831 ± 102 ± 44 

n *^ 10 + 0.020 

U.OiO_Q Q22 

1.9375-1.9500 

90± 11±2 

817 ± 104 ± 43 

0.298l°;°2? 

1.9500-1.9625 

80 ± 12 ± 3 

712 ± 105 ±42 

0.270t°°f2 

1.9625-1.9750 

91 ±12 ±2 

802 ± 105 ± 43 

0.2801^20 
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FIG. 21: The pp mass spectra in the mass regions near J/ij) 
(left) and ip{2S) (right). The curves are the the result of the 
fit described in the text. 


existing data with Mpp > 3 GeV/c^. Here A = 0.3 GeV 
and C is free fit parameter. It is seen that the asymptotic 
regime is reached at masses above 3 GeV/c^. 


VIII. THE J/tp AND ip{2S) DECAYS INTO pp 


The differential cross-section for ISR production of a 
narrow resonance (vector meson V), such as J/ip ^ decay¬ 
ing into the final state / can be calculated using [23] 


dcr(s, 9p) 
d cos 0* 


^ e+e-)B{V ^ f) 
my s 


W{s,xv,e;), 


( 12 ) 


where my and r(V —> e“'"e“) are the mass and elec¬ 
tronic width of the vector meson V , xy = 1 — royls, and 
B{V f) is the branching fraction of V into the final 
state /. Therefore, the measurement of the number of 
J/ip ^ pp decays in e+e“ —> ppy determines the prod¬ 
uct of the electronic width and the branching fraction: 
r(J/'0 ^ e'^e~)B{J/ip pp). The pp mass spectra for 
selected events in the J/ip and '0(25') mass regions are 
shown in Fig. 21. To determine the number of resonance 


events, both spectra are fitted with a sum of the prob¬ 
ability density function (PDF) for signal plus a linear 
background. The resonance PDF is a Breit-Wigner func¬ 
tion convolved with a double-Gaussian function describ¬ 
ing the detector resolution. The Breit-Wigner widths and 
masses for J/ip and 0(25) are fixed at the world-average 
values. The parameters of the resolution function are de¬ 
termined from simulation. To account for possible differ¬ 
ences in detector response between data and simulation, 
the simulated resolution function is modified by adding 
in quadrature an additional uq to both cr’s of the double- 
Gaussian function and introducing a shift of the central 
value of the resonance mass. The free parameters in the 
fit of J/0 mass region are the number of resonance events, 
the total number of nonresonant background events, the 
slope of background, ctg, and mass shift. In the 0(25) fit 
the UQ and mass shift values are fixed at those obtained 
for the J/0. 

The fit results are shown as curves in Fig. 21. Numer¬ 
ically, we find: Nj/.^ = 438 ±22 and fVy,(2S) = 22.2 ±5.7; 
the number of nonresonant events is 2'7 ± 8 for the 3- 
3.2 GeV/c^ mass interval and 7.9 ± 4.0 for the 3.6- 
3.8 GeV/c^ interval. These values are used to extract the 
nonresonant e“'"e“ ^ pp cross-section. Since the back¬ 
ground subtraction procedure for nonresonant events (see 
Sec. IV D) uses events with 30 < < 60, the mass spec¬ 

tra obtained with this cut may also be fit. The numbers 
of J/0 and nonresonant events are found to be 27 ± 6 
and 6 ± 4. The ratio of J/0 events with 30 < < 60 

to the number with Xp < 30, 0.061 ± 0.014 is in good 
agreement with value of /3pp-y = 0.048 ± 0.003 obtained 
in Sec. IV D. In the 0(25) mass region, no events are 
selected with 30 < yp < 60. The remaining fit param¬ 
eters are aa = 4.2 ± 1.8 MeV/c^ and Mj/.^ — = 

— (1.8 ± 0.7) GeV/c^. The fitted value of ac leads to a 
change in simulation resolution (11 MeV/c^) of 8%. 

The detection efficiency is estimated from MG sim¬ 
ulation. The event generator uses experimental data 
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for the angular distribution of protons in ij) ^ pp de¬ 
cays. This distribution is described by 1 -I- a cos^ -d with 
a = 0.660 ± 0.045 for J/V- [24, 25] and 0.67 ± 0.15 for 
'ip{2S) [28]. The model error in the detection efficiency 
due to the uncertainty of a is negligible. The efficiencies 
are found to be emc = 0.168±0.002 for J /'0 and emc 
= 0.161±0.003 for '0(25'). The data-MC simulation dif¬ 
ferences discussed earlier are used to correct the former 
efficiency values by (2.3±4.0)%. 

The cross-section for e+e” —> ipl VPl for 20° < 
0* < 160° is calculated as 

a(20° <ei< 160°) = 

yielding (11.0 ± 0.6 ± 0.5) fb and (0.57 ± 0.14 ± 0.03) fb 
for J/0 and 0(25), respectively. The radiative-correction 
factor R = cr/aBom, is 1.007±0.010 for J/0 and 1.011 ± 
0.010 for iIj{2S), obtained from a MC simulation at the 
generator level. 

The total integrated luminosity for the data sample is 
(232 ± 3) fb“^. From the measured cross-sections and 
Eq. 12, the following products are determined: 

r(J/0 ^ e+e-)B{J/'ip pp) = (12.0 ± 0.6 ± 0.5) eV, 



FIG. 22: The pp invariant mass spectrum for pp^ candidates 
in the (3.9-4.5) GeV/c^ mass range. The curves are the result 
of the fit described in the text. 


r(0(25) ^ e+e“)6(0(25) ^pp) = (0.70±0.17±0.03) eV, 

The systematic errors include the uncertainties of the 
detection efficiencies, the integrated luminosity, and the 
radiative corrections. 

Using the world-average values for the electronic 
widths [29], we calculate the 0 ^ pp branching fractions 
to be 

Bi J/'ip pp) = (2.22 ± 0.16) X 10“3 

and 

B{ili{2S) pp) = (3.3 ± 0.9) X 10“^ 

These values are in agreement with the corresponding 
world-average values: (2.17±0.08) x 10“^ [29] and (2.67± 
0.15) X 10-4 [29-31]. 


that r(y(4260) ^ e+e”) is much smaller than the cor¬ 
responding partial widths of all known = 1 char- 
monium resonances, while r(y(4260) —*■ J/07r“''7r-) is 
much larger [32]. The four-quark [33], hybrid [32, 34], 
meson- or baryon-molecular [35, 36] interpretations have 
been suggested to explain these unusual properties of the 
U(4260). Information about U(4260) decay modes other 
than J/07r+7r- can help clarify the nature of the Y (4260) 
resonance. In particular, charmless decays of the Y (4260) 
are expected in the hybrid model [32]. 

The pp mass spectrum is shown in Fig. 22 for ppy can¬ 
didates with pp mass in the range (3.9-4.5) GeV/c^. The 
mass spectrum is fit with the function 


dN dL ^ 
dm dm 


y/ao{m) -I- 


rnyFy 

ruy -I- irnyBy 


2 


+ B{m), 


(13) 


IX. UPPER LIMIT ON y(4260) ^ pp DECAY 

Recently, a resonant-like structure in the invari¬ 
ant mass spectrum of J/07r+7r- near 4.26 GeV/c^ 
was observed by BABAR in the ISR process e+e” —> 
J/tlnT~^'K~j [13]. This structure can be characterized 
by a single resonance with a width of about 90 MeV 
and is referred to as Y(4260). From the J/07r+7r- 
mass spectrum, the e+e” ^ Y(4260) ^ J/tpTr~^Tr~ 
cross section at the maximum of the y(4260) reso¬ 
nance was found to be (51 ± 12) pb. From the fact 
that the Y(4260) resonance is not observed in the to¬ 
tal e“''e- ^ hadrons cross section, one can conclude 


where cro(m) is the nonresonant cross section for ^ 
pp, B is the contribution of background processes, ay 
is the cross section at the maximum of the Y(4260) 
resonance, and my and Fy are the resonance mass 
and full width, respectively. The nonresonant cross 
section is described by Eq. 3 with the proton effec¬ 
tive form factor parametrized by the PQGD formula 
Fpp = C/m^ log^(m^/A^) with A = 0.3 GeV. The total 
background contribution in the (3.9-4.5) GeV/c^ mass 
range, mainly from e+e” —> pp7r°, is estimated to be 
7.5 ± 4.3 events. Due to the large uncertainty in the 
background, we cannot determine its slope from the fit; 
we assume a uniform background mass distribution. The 
mass and width of the Y(4260) are fixed at the values 
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obtained in Ref. [13]: uiy = (4259 ± 10) GeV/c^ and 
Fy = (88 ± 24) GeV/c^. The values of C, ay, and the 
interference phase (j) are free in the fit. The best-fit func¬ 
tion is shown in Fig. 22 as the dashed line. The solid line 
represents the fit with ay = 0 (null hypothesis). The 
optimal value of ay is 2.6 pb with a significance of 0.7 ct. 
The significance is estimated from the ratio of the val¬ 
ues of the likelihood function for the optimal fit and the 
fit to the null hypothesis. Since the best-fit value of ay 
is compatible with zero, we set an upper limit on the 
e+e“ —*■ y(4260) —> pp cross section. 

The mass spectrum is fit with different fixed values of 
the interference phase (j) and the upper limit at 90% confi¬ 
dence level (CL) is determined as a function of phase with 
the Neyman approach [37] using a Monte-Carlo tech¬ 
nique. The upper limit varies from 1.0 pb to 6.4 pb. The 
maximum value, corresponding to (j) = —7r/2, is chosen 
as a final result: ay < 6.4 pb at 90% CL. 

From the ratio of measured cross sections for 
y(4260) —*■ pp and F(4260) —*■ J/'k~ , we calculate 
an upper limit on the ratio of branching fractions: 


S(y(4260) ^ pp) 
6(y(4260) ^ J/i^TT+TT-) 


< 13% at 90% CL. 


(14) 


X. SUMMARY 

The process e+e" —> ppj is studied for pp invariant 
masses up to 4.5 GeV/c^. From the measured pp mass 
spectrum we extract the e“*"e“ —> pp cross-section and 
proton effective form factor. The form factor has a com¬ 
plex mass dependence. The near-threshold enhancement 
of the form factor observed in the PS170 experiment [10] 
is confirmed in this study. There are also two mass 
regions, near 2.25 GeV/c^ and 3 GeV/c^, that exhibit 
steep decreases in the form factor and cross-section. By 
analysing the proton angular distributions for Mpp be¬ 
tween threshold and 3 GeV/c^, the ratio |G_e/Gm| is ex¬ 
tracted. For masses up to 2.1 GeV/c^, this ratio is found 
to be significantly greater than unity, in disagreement 
with the PS170 measurement [10]. 

From the measured numbers of e“'"e“ ^ JPPl 
and e+e” —> ip(2S)j —> ppj events, the products 

r(J/V' ^ e+e-)B{J/'il) pp) = (12.0 ± 0.6 ± 0.5) eV, 


F(V’(25') ^ e+e-)B{'>p{2S) pp) = (0.70±0.17±0.03) eV, 


and their corresponding branching fractions are deter¬ 
mined: 


BiJ/ij} pp) = (2.22 ± 0.16) X 10“^ 


B{il){2S) pp) = (3.3 ± 0.9) X 10”^. 


The upper limit on Y (4260) —> pp decay is obtained at 
90% CL: 


B{Y ^ pp) 


< 13%. 


B{Y —> J/'i/lTT+TT ) 
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